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LONG-TERM  GOALS 

The  long-term  goal  of  this  project  is  the  development  and  application  of  new  methods  of  investigation 
for  the  use  of  Lagrangian  data  and  other  emerging  in-situ  and  remote  instruments  (drifters,  gliders,  HF 
radar  and  satellite)  that  provide  information  on  upper  ocean  advection.  Special  attention  is  given  to 
the  development  of  new  techniques  for  data  fusion  and  assimilation  of  data  in  Eulerian  numerical 
models. 

OBJECTIVES 

To  develop  methods  in  the  framework  of  the  ODD  AS  DRI  objectives,  namely  to  use  information  from 
drifting  sensors  to  improve  the  prediction  of  spreading  in  the  velocity  field.  The  project  has  the 
following  specific  objectives: 

1)  To  develop  new  techniques  for  assimilation  of  Lagrangian  data  and  of  other  tracers  in 
ocean  circulation  models. 

2)  To  develop  methods  to  improve  prediction  of  upper  ocean  advection  in  coastal  areas. 

3)  To  apply  and  test  the  methodologies  using  recently  collected  data  sets  in  the  framework  of 
Marine  Rapid  Environment  Assessment  (MREA)  experiments. 

APPROACH 

The  work  involves  a  combination  of  analytical,  numerical  and  data  processing  techniques.  The  method 
development  has  been  carried  out  in  collaboration  with  E.  Piterbarg  (USC),  A.  Molcard  (ESEET 
Universite’  of  Toulon),  A.  Poje  (CUNY),  V.  Taillandier  (CNRS,  Paris) 

WORK  COMPLETED 

1)  Pinal  revision  and  publication  of  a  paper  on  a  coastal  application  of  Lagrangian  data  assimilation 
using  drifter  data  (Taillandier  et  ah,  2008) 
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2)  Development  and  testing  of  a  method  for  traeer  data  assimilation  in  collaboration  with  L.  Piterbarg. 

3)  Participation  in  2007-2008  Marine  Rapid  Environmental  Assessment  (MREA)  experiments  in  a 
coastal  area  of  the  Mediterranean  Sea  (Ligurian  Sea,  Gulf  of  La  Spezia).  Indications  on  drifter  release 
strategy  have  been  provided. 

4)  Publication  of  a  paper  on  particle  spreading  based  on  MREA  results  using  drifters  and  HE  radar 
data  (Molcard  et  al.,  2008),  in  collaboration  with  A.  Molcard  and  P.Poulain. 

5)  Publication  of  a  paper  on  coastal  dynamics  in  presence  of  a  cape  (Magaldi  et  ah,  2008) 

RESULTS 

i)  Assimilation  and  velocity  estimation  from  Lagrangian  and  tracer  data. 

The  assimilation  methods  developed  during  the  present  and  previous  grant  period  for  Lagrangian  data 
are  based  on  correcting  the  model  velocity  field  by  minimizing  the  difference  between  the  observed 
positions  of  Lagrangian  instruments  and  the  model  forecasted  positions  from  numerical  trajectories.  In 
other  words,  the  Lagrangian  nature  of  the  data  is  taken  into  account  introducing  an  appropriate 
observational  operator  based  on  particle  advection. 

Applications  of  the  method,  based  on  a  variational  approach  that  optimally  blends  model  outputs  and 
Lagrangian  data,  have  been  performed  in  a  coastal  area  of  the  central  Adriatic  Sea  (Taillandier  et  ah, 
2008),  a  sub-basin  of  the  Mediterranean  Sea.  Data  from  surface  drifters  have  been  used  together  with 
outputs  from  a  state  of  the  art  model  (ROMS)  during  the  experiment  DOLCEVITA  in  2002.  The  area 
is  especially  challenging  since  it  is  characterized  by  a  swift  boundary  current  (Western  Adriatic 
Current,  WAC)  flowing  along  the  shelf  with  high  shear  and  by  an  interior  flow  with  high  variability 
and  intense  mesoscale  activity.  The  assimilation  has  been  applied  for  a  period  of  45  days  between 
October  1  and  November  15  2002,  when  there  is  good  drifter  coverage  of  the  WAC  in  the  vicinity  of 
the  Gargano  Cape,  providing  a  time  series  of  corrected  velocities  (estimates).  The  results  show  that  the 
assimilation  is  quite  effective.  The  error  on  the  particle  prediction  is  significantly  reduced  and  the 
estimates  of  residence  times  and  export  rates  from  the  WAC  appear  significantly  impacted  by  the 
analysis,  indicating  an  improvement  over  the  first  guess  with  respect  to  the  measured  quantities. 

An  approach  conceptually  similar  to  the  one  used  for  the  Lagrangian  data  assimilation,  has  been 
recently  developed  also  for  other  data  that  provide  advection  information,  such  as  tracer  data  from 
satellite  and  glider  data.  The  approach  takes  into  account  the  appropriate  observational  operators,  i.e. 
the  fact  that  the  tracer  evolution  can  be  described  as  a  forced  advection-diffusion  equation  and  glider 
positions  can  be  modeled  as  a  combination  of  forced  drift  and  advection.  Preliminary  work  on  the 
tracer  data  has  been  performed  in  collaboration  with  L.  Piterbarg.  The  method,  based  on  “confidence 
regions”  provides  velocity  estimates  by  combining  observation  of  a  continuously  distributed  tracer 
and  a  circulation  model  output,  and  is  presently  tested  using  the  twin  experiment  approach  with 
realistic  velocity  fields  from  a  Mediterranean  OPA  model.  The  area  of  interest  corresponds  to 
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Figure  1:  Example  of  velocity  correction  using  consecutive  information  from  tracer.  The  left  panels 
show  snapshots  of  velocity,  while  the  right  panels  show  the  tracer  distribution  after  three  days,  for 
the  real  (upper),  model  (middle)  and  corrected  (lower)  fields.  The  corrected  fields  are  significantly 

more  similar  to  the  real  ones  than  the  model  ones. 
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the  2007-2008  MREA  experiment  in  the  Ligurian  Sea,  where  in  situ  and  satellite  data  have  been 
eollected  and  are  available  for  future  applieations  and  real  data  testing. 

An  example  of  results  using  the  twin  experiments  is  shown  in  Fig.l.  A  tracer  is  released  from  a 
gaussian  initial  condition,  and  let  free  to  be  advected  by  the  "real"  velocity  field  for  three  days.  The 
tracer  distribution  is  observed  daily  and  the  model  velocity  field  is  corrected  using  consecutive  tracer 
information.  Snapshots  of  velocities  are  shown  in  Fig.l  (left  panels),  for  the  real  (upper),  model 
(middle)  and  corrected  (lower)  fields,  while  the  right  panels  show  the  tracer  distribution  after  three 
days.  The  corrected  fields  are  significantly  more  similar  to  the  real  ones  than  the  model  ones. 

The  method  can  be  envisioned  to  be  used  for  satellite  visible  or  SAR  data,  for  instance  in  the  case  of 
oil  spill  detection,  or  also  for  thermal  data,  provided  that  the  forcing  term  of  the  advection-diffusion 
equation  is  sufficiently  known.  We  are  presently  working  on  the  error  analysis,  starting  from 
theoretical  estimates  provided  by  Piterbarg  (2008),  and  improving  them  using  realistic  twin 
experiment  settings. 

ii)  Participation  in  MREA  experiments 

In  June  2007,  we  have  participated  in  a  Marine  Rapid  Environmental  Assessment  (MREA) 
experiment,  in  collaboration  with  NEIRC-NATO,  NRL,  OGS,  CNR,  INGV  (Italy)  and  LSEET- 
University  of  Toulone  (France).  The  experiment  called  LASIE  (Ligurian  Air  Sea  Interaction 
Experiment)  included  a  coastal  component  focused  on  the  Gulf  of  La  Spezia  (POET,  Predictive 
Oceanographic  Environmental  Trial).  A  further  MREA  experiment  is  presently  underway  (September- 
October  2008)  in  the  Ligurian  Sea,  funded  by  ONR  and  NASA.  The  experiment  is  mostly 
targeted  to  optical  measurements  in  the  frontal  area  of  the  Ligurian  current,  but  it  also  has  a  significant 
physical  oceanography  component.  In  particular,  a  glider  fleet  is  operating  during  the  experiment. 

During  the  experiments,  a  number  of  drifters  have  been  and  will  be  released  in  clusters  of  3-5  units, 
with  average  initial  distance  of  approximately  1  km  or  less.  We  have  participated  in  the  planning  and 
in  the  real  time  strategy  for  drifter  launching.  The  launching  in  tight  clusters  is  very  useful  for  the 
ODDAS-DRI  purposes,  since  it  allows  to  investigate  the  properties  of  drifting  arrays  and  of  their 
dispersion.  In  particular,  since  the  initial  distance  between  drifters  is  small,  it  provides  information  on 
processes  of  dispersion  at  various  scales,  from  submesoscale  to  mesoscale  and  up.  Two  main  areas 
have  been  chosen  for  successive  cluster  launches.  One  area  is  in  the  open  Ligurian  Sea,  next  to  a 
moored  buoy  (ODAS  buoy)  with  meteo-marine  sensors,  while  the  other  area  is  much  closer  to  the 
coast,  in  the  Gulf  of  La  Spezia.  In  this  way,  the  effects  of  coastal  dynamics  and  smaller  scale 
processes  can  be  specifically  investigated. 

We  are  also  presently  contributing  to  the  ongoing  MREA,  in  collaboration  with  NRL  providing 
diagnostics  based  on  the  NRL  NCOM  model  forecasts.  In  particular,  we  provide  maps  of  Einite  Size 
Layapunov  Exponents  (ESLE)  that  give  direct  information  on  regions  of  high  particle  spreading  and 
concentration,  therefore  indicating  fronts  and  barrier  structures.  These  information  can  be  of  interest 
for  glider  navigation,  signalling  possible  critical  areas  and/or  interesting  sampling  regions. 
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Figure  2:  Example  of  snapshot  of  VHF  radar  velocity  in  the  Gulf  of  La  Spezia  with  superimposed 
drifter  trajectories  (upper  panel).  Comparison  between  trajectories  of  drifters  (red  dots)  and 
synthetic  trajectories  computed  from  VHF  radar  (black  dots)  for  four  clusters  launched  in  the  Gulf 
In  all  cases  the  VHF  radar  based  trajectories  are  able  to  capture  the  cluster  spreading,  a)  and  b) 
correspond  to  launches  from  the  same  point  two  days  apart  (June  18  and  20)  showing  completely 

different  dispersion  behaviours 
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The  participation  to  the  MREA  experiments  has  been  very  fruitful  in  many  ways,  providing  us  with  a 
great  opportunity  to  test  methodologies  and  allowing  collaboration  with  a  number  of  investigators  and 
institutions.  They  include  NURC-NATO  (Michel  Rixen,  J.  Texeira,  Alberto  Alvarez),  NRL  (P.  Martin, 
R.  Allard,  G.  Peggion,  E.  Coelho),  OGS  (P.Poulain),  Elniversity  of  Toulon  (A.Molcard).  In  particular, 
the  collaboration  with  A.  Molcard,  is  sponsored  by  ONR  Global,  and  she  has  already  visited  RSMAS 
in  June  2007  to  work  on  data  analysis  and  planning  for  the  oncoming  experiments.  The  visit  was 
very  successful,  and  three  more  visits  are  planned  in  the  next  two  years.  Results  from  the  analysis  of 
the  MREA07  data,  performed  in  collaboration  with  A.  Molcard  are  presented  in  in). 

in)  Particle  spreading  and  analysis  of  drifters  and  HF  radar  data 

During  the  MREA07  coastal  component  POET  experiment,  a  VHE  Radar  WERA  data  set  (with 
resolution  of  250  m)  has  been  collected  in  the  period  15-30  June  by  A.  Molcard  (LSEET-University 
of  Toulone)  in  the  Gulf  of  La  Spezia  (Eig.2).  The  data  have  been  processed  and  analyzed  together  with 
coastal  surface  drifter  data  collected  by  P.Poulain  (OGS)  (Molcard  et  ah,  2008).  The  drifters  have 
been  launched  in  clusters  of  6-3  units,  aimed  at  investigating  the  significant  time  and  space  variability 
of  flow  dispersal.  In  particular,  two  clusters  launched  in  the  same  location  two  days  apart  (Pig. 2  a,b) 
show  strikingly  different  patterns  of  evolution.  The  first  cluster  is  swiftly  advected  with  little  particle 
spreading,  while  the  second  cluster  divides  into  two  branches  with  significant  drifter  spreading.  In- 
situ  drifter  trajectories  have  been  compared  with  synthetic  trajectories  computed  from  VHP  radar 
fields.  In  all  cases,  the  spreading  patterns  of  the  in-situ  clusters  are  correctly  reproduced  by  the 
synthetic  clusters,  including  the  dramatic  difference  between  the  two  launches  from  the  same  location 
(Pig,2).  The  positive  results  are  likely  to  be  due  to  a  combination  of  factors  and  primarily  to  the  high 
resolution  of  the  radar  fields,  and  they  confirm  that  radar  data  are  well  suited  for  the  study  of 
coastal  flows  in  limited  areas  with  complex  patterns  of  velocity  and  transport 

The  dynamical  reasons  for  the  high  variability  in  particle  dispersion  patterns  is  presently  under 
investigation,  considering  the  effects  of  local  and  large  scale  dynamics.  In  order  to  better  characterize 
the  transport  structures  inside  the  Gulf  we  computed  the  Pinite  Size  Lyapunov  Exponents  (FSLEs) 
based  on  the  radar  velocity  field.  FSLEs  are  a  diagnostics  indicating  areas  of  highest  dispersion  (blue- 
turquoise  lines  in  Fig. 3)  and  concentration  (red-yellow  lines).  These  lines,  or  ridges,  separate  different 
Lagrangian  structures  and  are  expected  to  act  as  particle  barriers.  During  the  first  launch,  occurred  in 
June  18  and  characterized  by  a  tight  cluster  evolution  (fig.  2a),  the  FSLE  field  is  characterized  by 
dominant  red-yellow  lines  (Fig. 3,  upper  panel)  that  propagate  and  separate  the  Gulf  interior  from 
the  open  sea.  Visualization  of  the  time  evolution  of  drifter  positions  superimposed  on  the  FSLE  fields 
shows  that  the  drifters  indeed  do  not  cross  the  red-yellow  lines  and  exit  the  Gulf  flowing  along  them. 
The  situation  is  very  different  during  the  second  launch  on  June  20,  characterized  by  high  cluster 
spreading  (Fig.2b).  In  this  case,  the  FSLE  field  (Fig. 3  lower  panel)  does  not  show  major  ridges  or 
barriers.  A  hyperbolic  point  can  be  detected  inside  the  Gulf,  indicated  by  the  crossing  of  blue  and  red 
lines.  Indeed,  the  drifter  trajectories  appear  to  spread  at  that  point,  following  the  out-flowing  red  lines 
in  the  two  southward  and  northward  directions. 

In  order  to  investigate  possible  connections  with  the  large  scale  dynamics,  we  considered  results  from 
NCOM  model  simulations  in  the  Ligurian  Sea  performed  by  G.  Peggion  during  the  MREA07 
experiment.  The  first  launch  of  June  18  occurred  during  a  south-eastern  wind  (Scirocco)  episode,  that 
reinforces  the  cyclonic  boundary  current,  as  shown  in  Fig. 4  (upper  panel).  This  is  expected  to  impact 
the  circulation  inside  the  Gulf,  since  the  boundary  current  partially  penetrates  the  Gulf  as  a  strong  and 
coherent  jet,  isolating  the  inner  part.  During  the  second  launch  on  June  20,  the  wind  was  weak. 
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Figure  3:  Examples  of  FSLE  fields  based  on  VHF  radar  data  in  the  Gulf  of  La  Spezia.  Red-yellow 
(blue-turquoise)  lines  indicate  regions  of  high  concentration  (dispersion).  Upper  (lower)  panel 
corresponds  to  June  18  (20),  when  the  tight  (high  spread)  cluster  in  Fig.  2a  (2b)  was  launched.  A 
dominant  ridge  or  barrier  (red-yellow  line)  is  present  in  the  upper  panel  at  the  entrance  of  the  Gulf, 
while  in  the  lower  panel  the  main  feature  is  a  hyperbolic  point  inside  the  Gulf  determined  by  the 

intersection  of  blue  and  red  lines. 
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Figure  4:  Velocity  fields  in  the  Ligurian  Sea  during  the  MREA07  experiment  from  NRL  NCOM 
model  (courtesy  of  G.  Peggion).  The  upper  panel  shows  the  response  to  Scirocco  (South-eastern 
wind)  wind,  with  intensified  boundary  current,  while  the  lower  panel  shows  the  response  to  Libeccio 
(South-western)  wind,  with  weaker  and  less  organized  boundary  current. 
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transitioning  from  the  Scirocco  event  of  June  18-19  to  a  Libeccio  (South  western)  event  oecurring 
after  June  22  (Fig.4,  lower  panel).  During  the  transition,  the  Ligurian  boundary  eurrent  is  less  strong 
and  organized,  and  this  eondition  is  likely  to  lead  to  a  less  eoherent  eireulation  also  inside  the  Gulf, 
that  becomes  more  directly  influenced  by  local  dynamics.  The  large  scale  dynamics  and  its  impact  on 
coastal  processes  are  under  further  investigation,  using  also  the  MREA08  results  where  FSLEs  are 
presently  computed  in  real  time. 

IMPACT/APPLICATIONS 

The  results  on  Eagrangian  data  assimilation  have  a  significant  impact  for  operational  systems  and  they 
provide  additional  value  to  drifters  and  gliders  especially  for  coastal  applications.  The  results  with 
FSLEs  based  on  VHP  radar  velocities  indicate  that  the  FSLE  technique  is  valid  also  in  presence  of 
small,  highly  variable  structures. 
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